Branching of quasinormal modes for nearly extremal Kerr black holes 
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We show that nearly extremal Kerr black holes have two distinct sets of quasinormal modes, which 
we call zero-damping modes (ZDMs) and damped modes (DMs). The ZDMs exist for all harmonic 
indices I and m > 0, and their frequencies cluster onto the real axis in the extremal limit. The DMs 
have nonzero damping for all black hole spins; they exist for all counterrotating modes (m < 0) 
and for corotating modes with < /i < /^c = 0.74 (in the eikonal limit), where jj, = m/{l + 1/2). 
When the two families coexist, ZDMs and DMs merge to form a single set of quasinormal modes as 
the black hole spin decreases. Using the effective potential for perturbations of the Kerr spacetime, 
we give intuitive explanations for the absence of DMs in certain areas of the spectrum and for the 
branching of the spectrum into ZDMs and DMs at large spins. 



PACS numbers: 04.25.-g, 4.25.Nx, 04.70.Bw 

I. Introduction. Nearly extremal Kerr (NEK) black 
holes (BHs) — i.e., BHs for which the dimensionless angu- 
lar momentum a 1 in the geometrical units, G = c = 
M = 1, used in this paper — have drawn much attention 
recently. Besides the mounting evidence for fast-rotating 
BHs in astronomy [T], NEK BHs have considerable the- 
oretical significance, e.g., in studies of weak cosmic cen- 
sorship [2 and in calculations of black-hole entropy [3] . 

For extremal Kerr BHs (a = 1) the near-horizon ge- 
ometry reduces to AdS2X jl]. This observation led 
to the Kerr/CFT conjecture, which states that extremal 
Kerr BHs are dual to the chiral limit of a two-dimensional 
conformal field theory [5]. In the past few years the ex- 
tremal Kerr spacetime and spacetimes violating the Kerr 
bound were shown to be unstable [Sj. The stability of 
BHs depends on the sign of the imaginary part of their 
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FIG. 1: (Color online.) Phase diagram for the separation 
between the single- and double-branch regime for NEK BHs. 
Large purple dots and gold crosses correspond to (Z, m) pairs 
with only ZDMs for perturbations with spin —2 and 0, re- 
spectively. Smaller blue dots correspond to {l,m) pairs with 
both ZDMs and DMs. The green line is the phase boundary, 
computed using the eikonal approximation. 



complex free vibration modes, called quasinormal modes 
(QNMs) [7j. Therefore the NEK QNM frequencies stud- 
ied here can shed light on the onset of extremal Kerr in- 
stabilities and prove useful in quantum field theory (for 
example, in the calculation of two-point functions 8]). 

Detweiler first used an approximation to the radial 
Teukolsky equation for NEK BHs (see also iOJ) to show 
that QNMs with angular indices I — m have a long decay 
time [1^. Using Detweiler's result, Sasaki and Nakamura 
PTj calculated QNM frequencies analytically and Ander- 
sson and Glampedakis proposed long-lived emission from 
NEK BHs [12 . However, there remains a long-standing 
controversy in the literature about what set of QNMs de- 
cay slowly [13) . whether long-lived radiation is possible 
[14] . and whether the imaginary part of the QNM fre- 
quencies vanishes as a — >■ 1 (compare [TT1[T3] with [13] )■ 
Despite the importance of this problem, our present un- 
derstanding of the QNM spectrum of NEK BHs is incon- 
clusive. 

In a recent paper jTHj , some of us used a WKB analysis 
to relate Kerr QNMs in the eikonal limit to spherical 
photon orbits around Kerr BHs. We pointed out that a 
subset of spherical photon orbits of extremal Kerr BHs 
reside on the horizon and that the corresponding QNMs 
have zero damping. This happens when the parameter 
^l = m/{l + 1/2) > Mc - 0.74. Hod [H] computed Hc 
in the eikonal limit, finding an approximate analytical 
result in agreement with [15] . 

In this work, we will show that the NEK geometry has 
two distinct sets of QNMs: zero-damping modes (ZDMs) 
and damped modes (DMs). ZDMs are associated with 
the near-horizon geometry of the BH, and they exist for 
all allowed values of / and m > (we classify modes using 
Leaver's conventions [T7], but we use units in which the 
BH has mass M = 1). DMs are associated with peaks 
of the potential barrier; in the eikonal limit, they exist 
when ^ < 0.74. This implies that ZDMs and DMs co- 
exist if < /I < 0.74. Figure [T] is a "phase diagram" in 
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QNM space, showing the regions where either the ZDMs 
or both the DMs and the ZDMs exist for scalar and 
gravitational perturbations with I < 15. We will dis- 
cuss this phase diagram further below. When the ZDMs 
and DMs coexist, and when the BH spin a is small, for 
each {I, m) there is only a single set of QNMs charac- 
terized by the overtone number n (where modes with 
larger n have stronger damping). For larger a, this set 
of QNMs appears to break into two branches. The DM 
branch originates from lower-overtone modes at smaller 
a, and its modes retain a finite decay rate as a — > 1. 
The ZDM branch originates from higher-overtone modes 
whose imaginary part becomes smaller than that of DMs 
as a — > 1, thereby forming the second branch. This is 
similar to the behavior of eigenmodes in quantum me- 
chanics when we parametrically split a single potential 
well into two potential wells (cf. Fig. [2] below, as well as 
[l8] for a somewhat analogous phenomenon in the theory 
of oscillations of ultracompact stars). 

II. Matched expansions. For e = 1 — a <C 1 and 

Lo — to/2 ^ 1, the radial Teukolsky equation can be writ- 
ten in a self-similar form when (r — 1) <C 1 and in an 
asymptotic form (by setting a = 1) when (r — 1) 3> ^/e 
(cf. [ini Uni). The solutions of the Teukolsky equa- 
tion in these regions (hypergeometric and confluent hy- 
pergeometric functions, respectively) can be matched at 
^Je {r — \) <^ 1 to provide the following condition for 
QNM frequencies: 

-7^S-2;S\nM^iS\n{8e) T^{2i5)T{l/2 + S - im ~ id) 

r2(-2i(5)r(l/2 + s- im + i5) 
r(l/2 - s-im- iS)r[l/2 + i{m-5- V^Cb)] 



r(l/2 -s-im + i5)r[l/2 + i{m + S- V^Cb)] 



(1) 



Here we denote the eigenvalues of the angular Teukol- 
sky equation by s^/m, and we define 6^ = Ira^ /A — (s + 
1/2)2 _ ^^j^ and Co = {lj — mQ.H)/\/t [note that VIh — 
a/ {ri^+a?) is the horizon frequency and r+ = — a? 

is the horizon radius]. Scalar, electromagnetic, and grav- 
itational perturbations correspond to spin s = 0, —1, —2, 
respectively. If we choose the conventions that TZe{5) > 
and Im{5) > when S'^ is positive and negative, respec- 
tively, then the left-hand side of Eq. ([T]) is usually a very 
small number, except when it is near the poles of the F- 
functions in the numerator. When m > 0, we can always 
find the solution near the poles at negative integers: 



1/2 + i{m - (5 - \/2aj) w -n, 
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FIG. 2: (Color online.) Plot of the potential term for a = 1, 
Eq. for p. = 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 (black-solid, red- 
dashed, red-dotted, blue-dotted, blue-dashed, and magenta- 
solid curves, respectively). The transition from single-branch 
to double-branch happens between ^ = 0.7 and [i, = 0.8. 



index of Kerr QNMs at lower spins. This set of solutions 
was first discovered by Hod [13] . The matched-expansion 
derivation shows that this set of modes depends on the 
near-horizon region of the Kerr BH. Equation (|3]) is quite 
accurate when |(5| 3> 1, but when \8\ < 1 it needs an ad- 
ditional correction [19]. However, the -y/e scaling of the 
decay rate is still correct when |^| < 1. The solutions to 
Eq. ([3| with to < are those that arise from the sym- 
metry w/_„i = — there are no solutions with to < 
and TZe{uj) > 0, when uj — m/2 is not small. Thus, the 
ZDMs only exist in the corotating regime to > 0. 

Another set of solutions of Eq. ([ij may exist when 
(5^ < and 2i5 sa — n, with n a positive integer. A more 
detailed analysis shows that, in this case, two nearly de- 
generate hypergeometric functions have comparable con- 
tribution to the near- horizon solution [19]. As a result, 
Eq. ^ is no longer valid when 2i5 ~ —n. As a consis- 
tency check, we looked for solutions with 2i5 ~ —n using 
Leaver's method and we did not find any. 

III. WKB analysis. The matched-expansion method 
assumes that uj w to/2, but Eq. ([T]) does not hold for 
modes which do not meet this requirement (i.e., DMs). 
To compute these modes, we will instead use a WKB 
analysis in the eikonal limit ^ ^ 1. The radial Teukolsky 
equation when Z 1 is [15] 
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Note that the overtone index n of these ZDM frequen- 
cies need not correspond precisely to the same overtone 



K ^ ~uj{r^ + a'^) + am, — ^ , 

dr* + dr 

\\[^ ^ Ai^ + a^u? - 2aTOa;, A = - 2r + . (4b) 
We define cj = w/j — iw/, and we note that the real and 



3 



imaginary parts scale as ujji cx I and w/ oc f , while the 
angular constant scales as Aim oc P. We only keep the 
leading-order terms in the eikonal limit in the following 
discussion (therefore all s-dependent terms are neglected, 
and the Aim are real). In Fig.[2]and below, we will refer to 
—Vr as "the potential". According to the WKB analysis 
and its geometric correspondence in |15j . the position 
of the peak of the potential asymptotes the horizon as 
a — >■ 1 for some of the corotating modes. For this set 
of QNMs, one can verify that V^' (where primes denote 
derivatives with respect to r*) scales as A^; thus, the 
peak r-Q of the potential becomes broad as ro approaches 
the horizon. It then follows that w/ oc yJV^.' / d^jVr — 0, 
and lor — m/2 in order to satisfy Vr{ujR,rQ) = for 
these modes. Assuming that tq = 1 + cy^ for the nearly 
extremal modes, where c is some constant, we can apply 
the eikonal equations in [15! and obtain 
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with J-Q = ■\/7m^/4 — Aimi'-o = m/2). Comparing this 
result with Eqs. ([S]) and ([5|, we can see the two sets 
of frequencies are essentially the same modes, although 
obtained in very different ways. Here J^q and 5^ differ 
by 1/4, which is reasonable because in the eikonal limit 
J-Q oc I and 5 oc I (making 1/4 a higher-order correction) . 

To build intuition about and 5, we look at Vr for 
extreme Kerr BHs, with lo replaced by m/2: 



Vr = 



(r2 + If 



(6) 



where L = ^-1-1/2 and a{ii) = Aim/L'^- According to the 
WKB analysis of the radial Teukolsky equation [2nj, the 
QNM frequencies are determined by the peak of the po- 
tential. As shown in Fig. [2] when ^ is large the maximum 
of the potential is at the horizon, r = 1, as expected for 
ZDMs. As /i decreases and falls below some critical value 
/ic, the peak moves outside the horizon, and the horizon 
becomes a local minimum of the potential. At the peak 
LUi is nonzero because (fiVr / dr1\ra ^ 0, so we have DMs. 
The criterion for having no peak outside the horizon is 
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fi^ — a{fi) + -fi^ > for r = 1 



(7) 
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i.e., J"o > (or > 0). The values at which J^q (or S"^) 
vanish lead to the condition for the critical /ici a{fic) = 
jHc- If we use the approximation a(/i) « 1 — a^uj'^[l — 
/i2)/(2L2) [TS], this will reproduce Hod's approximate 
analytical result fic ~ [(15 - 7l93)/2]i/2 [jg. We can 
obtain the exact fic (in the eikonal limit) by inserting 
a{fic) = into the Bohr-Sommerfeld condition for a 
derived in ilS : 
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+ ^cos2 0=(l-|/.|)7r, (8) 



where 9+ = it— 9- and 0_ = arcsin(-\/3— 1) are the angles 
at which the integrand vanishes. Therefore we have 
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which yields the numerical value /ic — 0.74398. In the 
eikonal limit, when ^ > NEK BHs have only ZDMs 
("single-phase regime") ; when < fi < fic, both DMs 
and ZDMs exist ("double-phase regime"). 

IV. Phase boundary. Although there is a clear cri- 
terion for determining the boundary between the single- 
phase regime and the double-phase regime in the eikonal 
limit (when fi < fj.c, the peak of the potential no longer 
resides on the horizon) it is not immediately clear if a 
similar criterion holds when I is small. For scalar per- 
turbations, however, we can write the radial Teukolsky 
potential for extreme-Kerr BHs with generic Z,m, under 
the assumption that w = to/2 (and, therefore, the oAim 
remain real for the ZDMs): 
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It is not difficult to see that there is still no peak out- 
side the horizon when |to2 > oAim, or 7"q > 0. For 
electromagnetic and gravitational perturbations the po- 
tential terms —Vr are complex functions, thereby mak- 
ing the positions of their extrema more difficult to define. 
Detweiler [2T], however, has shown that the radial func- 
tion can be transformed so that it satisfies a differential 
equation with a real potential. Using this potential, the 
criterion to exclude peaks outside the horizon is |19j : 



J^t = Jto- 



uj=—j>0. (11) 



Note that this expression respects the pairing symmetry 
-sAim — sAim + 2s, and that for all s, and 5^ differ 
from each other only by 1/4. For s = 0, —2 and 2 <l < 
100, we have searched all QNMs numerically and have 
not found any mode simultaneously satisfying 5^ < Q and 
J-'2 > 0; therefore, the sign of J2 also determines whether 
a peak exists outside the horizon. In addition, we have 
used Leaver's continued-fraction algorithm to determine 
the phase boundary numerically. As shown in Fig. [ij the 
actual phase boundary matches the criterion predicted 
by the eikonal limit, fj, — fi^- In addition, for scalar 
and gravitational perturbations, we find numerically that 
modes are in the single-phase regime when J-"2 > for all 
I < 15. This reinforces our physical understanding that 
DMs are associated with a peak of the potential outside 
the horizon, while ZDMs are somewhat similar in nature 
to the s-modes in ultracompact stars |18j . 
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FIG. 3: (Color online.) QNM frequencies with I = 10 for 
NEK BHs. Contours are constant values of the logarithm 
of the continued fraction in the complex plane; darker shad- 
ing indicates values near zero. The + symbols are the ZDM 
predictions, a x is the lowest-overtone WKB prediction from 
jl5| . and the box is centered at the WKB prediction from 
Eq. ([5|. No branching is observed for modes with m — 8. 
Note that the closed contours with light shading have large 
values and do not correspond to any QNM. Further discussion 
of the figure is in the text. 



V. Bifurcation. Schwarzschild and slowly spinning 
Kerr BHs have a single set of QNMs for each I, m that are 
characterized by their overtone number n. If the ZDMs 
originate from modes at higher-n than the DMs when the 
BH spin is low, then when the spin increases beyond a 
critical value = 1 — eg, a, single set of QNMs may split 
into two branches. 

We numerically investigate this bifurcation effect by 
examining the complex QNM frequency plane to search 
for solutions of Leaver's continued-fraction equations 
[T71 [TO] . In Fig. [3] we plot the contours of constant 
value of the logarithm of the continued-fraction expan- 
sion, truncating at = 800 terms. The QNM frequen- 
cies correspond to the local minima of this sum, where 
the contours cluster. The shading indicates the value of 
the fraction, with darker values nearly zero. 

When /X < /Xc, a single set of QNMs splits into two 
branches for increasing a (see the left-hand panels of 
Fig. |3j where Z = 10, m = 7, as the spin increases from 
a = 0.9990 to a = 0.9999 from the upper panel to the 
lower). The ZDM branch is quite accurately described 
by Eq. ([3|; the imaginary part of the ZDMs scales like 
^/e, and they move towards the real axis as e ^ 0. The 
DM branch changes relatively little with increasing spin 
(it is expected that the WKB peak can only support a 



finite number of modes [T5], and there are only 3 DMs 
in the lower- left panel). In this case, the WKB formulae 
of [15] are in good agreement with the lowest-overtone 
DM (marked with a x in the figure). 

For fi > fic there is no bifurcation, and the modes are 
predicted fairly well by Eq. ([s]). We can see this in the 
right-hand panels of Fig. [3] where I — 10, m ^ 8 and we 
again raise the spin from a = 0.9990 to a = 0.9999. For 
the TO = 8 modes, we also mark the leading-order WKB 
prediction of Eq. ([s]) with a box. For the bifurcation 
effect, we can define a benchmark Oc = 1 — ec as the 
BH spin at which the imaginary part of the fundamental 
ZDM is equal to that of the fundamental DM: 



2V2 



(1 + 2|<5|) = 



2 d^Vr 



(12) 



The right-hand side of Eq. ( 12 ) can be evaluated using 



the approximate WKB formula in [T3]. Since both sides 
of Eq. ( 12 ) depend on e, we solve for iteratively; this 
converges quickly for a variety of initial spins. By com- 
puting Ec for I < 15 and < to < (/ + l/2)/ic, we find 
that L'^ec = 10-^(11.6 - 3.12^ - 18.0^^) is a reasonable 



fitting formula. For the I — 10, m — 7 case, Eq. ( 12 ) 
gives ec ^ 10~^, which is in agreement with numerical 
results; for the / = 2, to = 1 case it gives ec ~ 10"'^. 

In Fig. [3] however, it is clear that the bifurcation actu- 
ally starts when the fundamental ZDM's imaginary part 
equals the imaginary part of the highest-overtone DM (in 
Fig. |3]it is the third overtone). This happens at a spin 
fls < Uc- Because we do not have a good estimate of 
the number of modes in the DM branch (beyond the fact 
that it should be proportional to L and a function of fi 
in the eikonal limit [19]) and because WKB techniques 
are not accurate for these high-overtone DMs, finding an 
analytic solution for a<, remains an open problem. 

VI. Conclusions. We identified two different regimes in 
the NEK QNM spectrum. In the double-phase regime, 
we found that the lowest ZDM becomes less damped than 
the lowest DM at some critical Uc, for which we provided 
an analytical estimate. For sufficiently large a, Eq. ^ is 
accurate at the least for those ZDMs with smaller decays 
than the point where the branches bifurcate. We estimate 
that the number of ZDMs below the bifurcation is cx 
\/es/e [in]- In the future, we would like to investigate 
the behavior of the ZDM branch in the high-overtone 
limit |22j . where these approximations break down. 
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